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Using the post-Newtonian (PN) expansion technique of the gravitational wave pertur¬ 
bation around a Schwarzschild black hole, we calculate analytically the energy flux of grav¬ 
itational waves induced by a particle in circular orbits up to the 5.5PN order, i.e. 0{v^^) 
beyond Newtonian. By comparing the formula with numerical data, we find that the error 
of the 5.5PN formula is about 4% when the particle is on the last stable circular orbit. We 
also estimate the error AN in the total cycle of gravitational waves from coalescing compact 
binaries in a laser interferometer’s band produced by using the post-Newtonian approxima¬ 
tions. We hnd that, as for the neutron star-black hole binaries, the 4.5PN approximation 
gives AN < 1 for a black hole of mass M < 4OM0, while it gives AN > 1 for a black hole 
of mass M > 40M©. 


§1. Introduction 

Gravitational waves from coalescing compact binaries are the most promising 
candidates which will be able to be detected by the near-future, ground based laser 
interferometric gravitational wave detectors such as LIGO, VIRGO, TAMA, GEO600 
etc.0^ 

If a neutron star or a small black hole spirals into a black hole with mass 
< 300Mq, the inspiral wave form will be detected by the above detectors. When 
a signal of gravitational waves is detected, we will try to extract parameters for 
binaries, such as masses and spins etc., from inspiral wave forms using the matched 
filtering technique. 13 In this method, parameters for binaries are determined by 
cross-correlating the noisy signal from detectors with theoretical templates. If the 
signal and the templates lose phase with each other by one cycle over 10^ - 10^ 
cycles as the waves sweep through the interferometer band, their cross correlation 
will be significantly reduced. This means that, in order to extract the information 
optimally, we need to make theoretical templates which are accurate to better than 
one cycle during an entire sweep through the interferometer’s band. 13 

The standard method to calculate inspiraling wave forms from coalescing binaries 
is the post-Newtonian expansion of the Einstein equations, in which the orbital 
velocity v of binaries is assumed to be small compared to the speed of light. Although 
the post-Newtonian calculation technique is being developed to apply to higher order 
calculation,® it becomes more and more difficult and complicated. Thus, it would 
be very helpful if we could have another reliable method to calculate the high order 
post-Newtonian corrections. 

As an alternative method, the post-Newtonian expansion of the black hole per- 
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turbation has been developed. There, one considers gravitational waves from a 
particle of mass /r orbiting a black hole of mass M, assuming ^ M. Although 
this method is restricted to the case /r <C M, we can calculate very high-order 
post-Newtonian corrections to gravitational waves by means of a relatively simple 
analysis, compared to the standard post-Newtonian analysis. 

Since LIGO and VIRGO will be able to detect gravitational waves from binaries 
of mass less than ~ SOOMq, it is important to construct templates for such binaries. 
The frequency of gravitational waves for such massive binaries, however, comes into 
the frequency band for LIGO and VIRGO at r/M ~ 16(1OOM0/M)^/^, i.e., the 
highly relativistic region. We do not know whether the convergence property of the 
post-Newtonian approximation is go(^ or not for such a highly relativistic orbit. As 
was reported by Tagoshi and Sasaki, the post-Newtonian convergence of the total 
orbital phase during the detectable frequency band for LIGO and VIRGO can be 
very slow for these binaries. Hence, it is an urgent problem to clarify the extent to 
which the convergence property of the post-Newtonian expansion is good. For this 
purpose, we study the energy loss rate of the binary up to order in this paper, 
where v := . 

This direction of research was first nursued done analytically by Poisson to 
O(u^) and numerically by Gutler et al.cP to A highly accurate numerical 

calculation was done by Tagoshi and Nakamura, 0) and they gave a fitting formula 
to 0(u®). They found the appearance of Inu terms in the energy flux at 0(u®) and 
at 0(u®). They also clarified that an accuracy of the energy flux to at least 0(u®) 
is needed to construct template wave forms for coalescing binary neutron stars. An 
analytical calculation to the same oiHer was done by Tagoshi and Sasaki 0*, based on 
the formulation built up by Sasaki, 0^ and they confirmed the result of Tagoshi and 
Nakamura. Recently, a systematic method to calculate the higher order corrections 
was developed by Mano, Suzuki and Takasugi.l^ 

Note that, both in these analyses and in this paper, the particle is treated as 
a test particle. According to the order of magnitude estimation, the quadrupole 
moment of a star produced by the tidal deformation of the star affects gravitational 
waves at O(u^°).0 Thus we cannot regard the realistic compact star as a test particle. 
Nevertheless, we ignore this fact and treat the star as an ideal structureless particle 
in this paper, expecting that we can still obtain useful information by analyzing such 
an ideal system. 

This paper is organized as follows. In §2, we give the general formulas and 
conventions used in this paper very briefly. The full description of the formalism 
can be found in Tagoshi and Sasaki® and references therein. In §3, we show the 
expressions of the energy flux extended to 0{v^^). In §4, we discuss the accuracy of 
the post-Newtonian expansion. §5 is devoted to a summary. 

Throughout this paper we use the units such that c = G = 1. 


§2. General formulation 


For the reader’s convenience, we attach a minimum amount of the explanation 
about the black hole perturbation formalism in order to make the convention defi- 
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nite. For the reader who would like to know the details, see the reference(Some 
typographical errors in the formulas in this reference are corrected here). We con¬ 
sider the case in which a particle of small mass ^ travels a circular orbit around a 
Schwarzschild black hole of mass M ^ fi. 

To calculate the gravitational luminosity, we consider the inhomogeneous Teukol- 
sky equation. 


Z\^ 


iA. 

dr 


AA. 

A dr 


+ U{r) 


PurtiLd (^) — Tlmuj (j ")) 


( 2 - 1 ) 


where 


^2 ^ 


U{r) = —^Lo^r^-Aiuj{r-3M)^ -{£-!){£+ 2), A = r{r-2M), (2-2) 

and is the source term which reflects the energy momentum tensor of the small 
particle. We omit the explicit form of T(^muj{r) here. 


We solve Eq. (2T) using the Green function method. For this purpose, we need 
a homogeneous solution of Eq. (^) which satisfies the boundary conditions 


Dii^A'^e for r* —> — oo, 

+ r~^ for r* -|-oo. 


(2-3) 


where r* = r + 2Mln(r/2M — 1). Then the outgoing-wave solution of Eq. (O at 
infinity with appropriate boundary conditions at horizon is given by 


Rimu! ir 


„6„iLor poo 

/ drRZTe^.{r)A-^ 

2lU)BZ J2M 

Am.. 


(2-4) 


In the case of a circular orbit, the specific energy E and angular momentum L of 
the particle are given by 


E = {ro- 2M)/\JrQ{rQ - 3M), (2-5) 

and 

L = AMro/ - 3M/ro, (2-6) 

where ro is the orbital radius. The angular frequency is given by 17 = 

Dehning shim by 

ohm = h)£{£ + 1)(^ + ZGim o) Ero/(ro - 2M), 

-ihm = [(^ -m + -lYim o) i/ro, 

-2&fm = -2T^m ^-^5 0^ Li7, (2-7) 

where nY£m{G,p) are the spin-weighted spherical harmonics, is found to take 

the form 


Zimui — 


( 2 - 8 ) 












4 


T. Tanaka, H. Tagoshi and M. Sasaki 


where 


^im — . 


TT 




-ohm - ‘^i-ibem ( 1 + ^wro/(ro - 2M)^ 


+i_ 2 bimi^ro{l - 2M/ro) ^ 1 - M/ro + -iwro 


1 


Z?' 

IL, 


'im 


+ 


i-ibim --2 bem (l + {ro - 2M)) ro-R^^Vo) 


+ 2-‘^^^fnroR}Z (ro) 


(2-9) 


In terms of the amplitudes Zim^ the gravitational wave luminosity is given by 

^=E E \Zem\V27ru;^ (2-10) 

® e=2m=-e 

where oj = mfi. We calculate this quantity in the post-Newtonian expansion, that 
is, in the expansion with respect to u = Thus the only remaining task 

is to calculate the series expansion of the ingoing-wave Teukolsky function in 
terms of e := 2Moj = O(u^) and z := ujr = 0(v) and the Wronskian in terms 
of e. 


§3. The gravitational wave luminosity and the accuracy of the 

post-Newtonian expansion 


As noted in the introduction, a new method to calculate the Teukolsky functicm 
in the post-Newtonian expansion was constructed by Mano, Suzuki and Takasugi.H^ 
However, we have derived the following results by using the formalism developed 
by Sasaki,™ and the details of derivation is too complicated and seem not worths 
of note here, because a new more powerful method has been developed. Thus we 
present only the final results of the gravitational wave luminosity. The formulas for 
B^^ and the luminosity of e^h mode are given in Appendix A. 

The total luminosity becomes astT 


dE 

dt 


dE\ 

dt J N 


. 1247 2 , 3 

1- V -b 47ru'^ 

336 


44711 4 

- V 

9072 


8191 TT . 
- 

672 


/6643739519 
y 69854400 


17127 1712 In u 16 vr^ 3424 In 2 \ g 

105 105 ~3 105 ) ^ 


16285 TT 7 

- V 

504 


/ 323105549467 2325977 232597 Inu 

I 3178375200 4410 4410 

1369 TT^ 39931 In 2 47385 In 3 

-b 


126 


294 


1568 


V 
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+ 


p65978667519 7 r 6848 7 vr 68481n?;7r 13696 7 r ln2^^ 9 

I 745113600 105 105 105 


2500861660823683 9166284677 916628467 Inu 424223vr^ 

+ 1 -+ 


2831932303200 


7858620 


7858620 


6804 


83217611 In 2 47385 In 3 

+ 


,10 


1122660 


196 


+ 


8399309750401TT 177293 7 vr 177293 In ?; 7 r 
+ —+ ■ 


101708006400 


1176 1176 

8521283 vr In 2 142155 vr In 3 


+ - 


.11 




N 


17640 784 

1 - 3.711309523809524 + 12.56637061435917 

-4.928461199294533?;'^ - 38.29283545469344?;^ 


+ (115.7317166756113 - 16.3047619047619 In ?;) 
-101.5095959597416?;^ + 


(-117.5043907226773 + 52.74308390022676 In ?;) 

+ (719.1283422334299 - 204.8916808741229 In ?;) 

+ (-1216.906991317042 + 116.6398765941094 In ?;) ?;^° 

+ (958.934970119567 + 473.6244781742307 In ?;) ?;^^ + • • • 


(3-1) 


Using the above results, we compare the formula for the gravitational wave flux 
with numerical results and investigate the accuracy of the post-Newtonian expansion. 

A high precision numerical calculation of gravitational waves from a particle in a 
circular orMt around a Schwarzschild black hole has been performed by Tagoshi and 
Nakamura. 0 Since no assumption was made about the velocity of the test particle, 
their results are correct relativistically in the limit /x <C M. In that work, dE/dt was 
calculated only for £ = 2 ~ 6. Then, for the orbital radius tq < lOOM, we calculate 
dE/dt again for all modes of £ = 2 ~ 6 and for i = 7 with odd m. The estimated 
accuracy of the calculation is about 10 “^^, which turns out to be accurate enough for 
the present purpose. As for the radius ro > lOOM, since it inexpensive to calculate 
dE/dt, we use the data calculated by Tagoshi and Nakamura™ which contain modes 
from £ = 2 to 6 . 

In Figs. 1 and 2, we show the error in the post-Newtonian formulas as a function 
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Figs. 1 and 2 The error of the post-Newtonian formulas as a function of the Schwarzschild coordinate 
radius r for 6 < r/M < 100 (Fig. 1) and 100 < r/M < 5000 (Fig. 2). In Fig. 1, contributions 
from £ = 2 to 7 are included. In Fig. 2, contributions only from £ = 2 to 6 are included in both the 
post-Newtonian formulas and in the numerical data. 

of the orbital radius r. We refer to the post-Newtonian formula which include terms 
up to the order u” as the (n/2)PN formula. The error of the post-Newtonian formula 
is defined as 



(3-2) 


where {dE/dt)n and (dE/dt) denote the (n/2)PN formula and the numerical result, 
respectively. In the plot of Fig. 2, only the contributions from £ = 2 to 6 are included 
in both the post-Newtonian formulas and the numerical data. We can see that, at 
small radius less than r ~ lOM, the error of the IPN and 2.5PN formulas are larger 
than the other formulas. On the other hand, the Newtonian and the 2PN formulas 
are very accurate within this radius. This is because those formulas coincide with 
the exact one accidentally at a radius between 6M and lOM. The error of each 
post-Newtonian formula at the inner most stable circular orbit, r = 6M, becomes 


as follows; 12% (Newtonian), 66% (IPN), 8.6% (1.5PN), 3.4% (2PN), 42% (2.5PN), 
11% (3PN), 5.4% (3.5PN), 17% (4PN), 8.4% (4.5PN), 6.5% (5PN), 4.1% (5.5PN). 


As is expected, the errors of the post-Newtonian formulas decrease almost linearly 
up to r ~ 5000M in a log-log plot. This fact also suggests that the numerical data 
have accuracy of at least 10 at ■p rsy 5000M. 

In order to examine exactly to what order the post-Newtonian formulas are 
needed to do accurate estimation of the parameters of a binary using data from laser 
interferometers, we must evaluate the systematic error produced by incorrect tem¬ 
plates, llil However, here we simply calculate the total cycle of gravitational waves 
from a coalescing binary in a laser interferometer band and evaluate the error pro¬ 
duced by the post-Newtonian formulas. It has been suggested that whether the error 
in the total cycle is less than unity or not gives a usrful guideliim to examine the 
accuracy of the post-Newtonian formulas as templates® (see alsoll®). 

We ignore the finite mass effect in the post-Newtonian formulas and interpret 
M as the total mass and // as the reduced mass of the system. Further, since the 
error of the 5.5PN formula for dE/dt is only a few percent even at r = 6M, and 
since the error decreases monotonically as the orbital radius increases, we regard the 
5.5PN formula as if it were the exact formula. 
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The total cycle N of gravitational waves from an inspiraling binary is given by 

17 dE/dv 




N= dv- , , , 

Jvf vr \dE/dt\ 


(3-3) 


where Vi = ijvl, Vf = and rj and rf are the initial and final orbital 

separation of the binary. In the test particle limit, which we assume, dE/dv is given 
by 

dE d [ 1 — \ 



As noted before, we regard the 5.5PN formula (3T) for dE/dt as if it were the exact 
formula. We then introduce P(u) = {dE/dt)/[dE/dt )n and(5(t') = {dE/dv)/{dE/dv )n ^ 
where {dE/dv)^ = —/uu, and define the approximated total cycle as 


= 


dv 


IVf 


17 {dE/dv)N Qn{v) 
TT {dE/dt)M Pn{v) ’ 


(3-5) 


where Qn{v) and Pn{v) are given by expanding Q{v) and P{v) in terms of v and 
truncating them at order v^. Note that we do not expand Qn{v)/Pn{v) in terms of 
V further in calculating the above integral because such an expansion will produce 
additional errors in The error in the total cycle by using the (n/2)PN formula 

is defined as 

^N^n) ^ jyin) _ 

Note that, because we have regarded the 5.5PN formula for dE/dt to be exact, 
is produced only by the error in dE/dv. 

The results are listed in Table I. We find that the post-Newtonian approximations 
with n < 5 produce errors much greater than 1. This slow convergence was observed 
in previous worksFor a system of a neutron star with mass 1.4M© 
and a black hole with mass much larger than 1.4M0, the convergence becomes much 
slower. Accidentally, the 3PN formula seems to be good irrespective of the masses 
of the binary. At n = 7 and 8, the post-Newtonian approximations produce AN ~ 1 
for a binary with mass of each star less than about IOMq but AN ~ 10 for the case 
where one of binary stars is a black hole with mass ^ 40 Mq. At n > 9, AN < 1 for 
all the cases investigated here. 

From these results, we can say that at least the 3PN formula is needed for tem¬ 
plate wave forms. The 4.5PN or higher order approximation seems to be sufficient 
for constructing template wave forms. It should be noted, however, if the black hole 
mass is greater than 70Mq, AN becomes greater than 2 even in the case of the 5.5PN 
formula. Nevertheless, this does not necessarily mean that the post-Newtonian ap¬ 
proximations, even if we include the 5.5PN terms, is inadequate for calculating the 
theoretical wave form for such a system. To obtain a definite conclusion, one must 
perform a_more detailed analysis of the systematic error as done by Cutler and 
Flanagan. tiJ) 

One reason for the relatively slow convergence of the post-Newtonian expansion 
in calculating the cycle N is that we are integrating to a small radius r = 6M, which 
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Table I. The error in the total cycle of gravitational waves from coalescing compact binaries 

in a laser interferometer’s band produced by using the (n/2)PN formula. The initial frequency 
is lOHz, which corresponds to the initial radius r/M ~ , and the final frequency 

is the one at which the binary is r = 6M. Note that is produced only by the error in 

the formula for dE/dv. 


o 

(1.4M0,1.4M0) 

(1OM0,1OM0) 

(1.4M0,1OM0) 

(1.4M0,4OM0) 

(1.4M0,7OM0) 

0 

-118 

12 

9.6 

100 

140 

2 

239 

67 

223 

312 

353 

3 

10 

6.3 

17 

40 

56 

4 

-1.3 

1.2 

1.6 

13 

24 

5 

11 

8.8 

22 

56 

78 

6 

-0.17 

-0.018 

-0.19 

0.75 

2.7 

7 

1.1 

1.1 

2.5 

7.9 

13 

8 

0.94 

0.91 

2.1 

6.7 

11 

9 

0.11 

0.12 

0.25 

0.93 

1.7 

10 

0.20 

0.20 

0.46 

1.5 

2.6 

11 

0.17 

0.17 

0.39 

1.3 

2.2 


implies a value of u ~ 1. Further, we can understand some of the above results by a 
simple order of magnitude estimation (see paper TSTS, section IV). Since the cycles 
are mainly accumulated around ~ lOHz, which is the lowest detectable frequency 
region of the laser interferometers such as LIGO and VIRGO, N is approximately 
given by 


V ~ 1.9 X 10' 


3 / lOMp y/a /M\ 

V M ) W)' 


(3-7) 


where M and /x are the total mass and reduced mass, respectively. This implies that 
the template must have an accuracy of at least 


~ 5 X 10“^ 


M 

lOM© 



(3-8) 


when the frequency of gravitational wave becomes lOHz. Then, for fixed M, we 
need much better accuracy in the case <C M than in the case of equal mass, since 
N becomes large. On the other hand, for fixed /x, the required accuracy seems to 
become lower as M becomes larger. However, since the non-dimensional radius r/M, 
at which the frequency of gravitational waves becomes lOHz, becomes smaller for 
larger M, the importance of higher order corrections increases. This is the reason 
why the convergence of the case (1.4M0,7OM0) is slower than that of (1.4M0,4OM0). 


§4. Summary 


We have calculated the energy flux of gravitational waves induced by a particle 
of small mass in circular orbits around a Schwarzschild black hole up to 0{v^^) 
beyond the quadrupole formula analytically. By comparing the analytically derived 
post-Newtonian formulas with numerical results, we discussed the accuracy of the 
post-Newtonian expansion. We have found that the accuracy of the 5.5PN formula 
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is good because the error of the dE/dt formula is about 4% even at the inner-most 
stable circular orbit r = 6M. 

We have also calculated the total cycle of gravitational waves in a laser interfer¬ 
ometer’s band and estimated the error produced by the post-Newtonian expansion. 
We have found that in the case when each star of a binary has the same mass, 
the 4.5PN approximation gives a sufficiently accurate cycle irrespective of the value 
of total mass. On the other hand, in the case of neutron star-black hole binaries 
for which the mass of the black hole greater than ~ 40Mq, the convergence of the 
post-Newtonian expansion is slow. When the total mass is about 40Mq, the error is 
about 1 for the 5.5PN approximation. When the total mass is greater than ~ 70Mq, 
the error becomes greater than 2 even if we go to the 5.5PN order. It is difficult to 
conclude, however, whether the post-Newtonian expansion can or cannot be used to 
construct theoretical templates for such systems only from the above results, because 
the value of AN is marginal. We need a detailed analysis of the systematic error, in 
the parameter estimation, produced by the post-Newtonian expansion, such as the 
one done in Ref. 12). This is left for future work. 
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Appendix A 

Formulas for Rff and and the gravitational wave luminosity 


The post-Newtonian expansion of the ingoing Teukolsky function in the near 
zone, where z is small, is given by 


R 


in _ 

2a; 


4 2:^ 8 i 

-^-^ 

5 15 


22 2i ^ 23z® , 2i 
105 ~ 35 ^ 1890 ^ ^ 


13 

41580 


i 59 z^^ i 83 

24948 ^ 12972960 2162160 ^ ~ 1945944000 


i 

277992000 



-8z3 3i 4 Sz® 13i « 109z^ 341i 

-z^-z“ H-^- 

5 5 63 90 1890 22680 


9403 z^ 293 i _ io 38963 z^^ 75529 i 12 ) , 

3118500 ~ 594000 ^ 567567000 9081072000 ^ P 


4z 2 ^ 123317 z^ , 231479 i ^^5 889954 z® 454499 i ^^7 

36750 110250 ^ 1157625 ~ 2315250 ^ 



10 


T. Tanaka, H. Tagoshi and M. Sasaki 


^ ;i5iU08iiz g 214 z’ 

5501034000 5501034000 ^ ^ 

1177 In z 107 i 7, 2461 z® In 

- 1 - 2 : In 2 ;- 

11025 3675 396900 


11025 


yy220 


I -66823 z^ 99851 i 4 504569 2 ® 2488639 ? 

’ I 12250 55125 ^ 694575 “ 3969000 

107 ? 4 , 4282 ® In 2 1391 ? e , \ o 

H- 2 In 2 H-^- 2 '’ In 2 e'* 

350 6615 18900 j 

U71A87z‘^ 263 ? 3 214 2 ^ ln 2 \ 4 

I 220500 ~ TM ^ 525 j ^ ’ 


99851 ? 4 


2488639 ? g , 428 2 ® In 2 
3969000 ^ 525 


/ 471487 2 ^ 
y 220500 

^42® 2 ? 0 

- 


2^ In 2 ' 


? 3 29 2 ® ? 

20790 4620 


■ 10 472 ®^ 

1621620 


? 42 I 23 2 ®® ^ ? 

'294840 ^ 64864800 29937600 ' 


-10 2 ^ 53 ? . 2 ® 

- — - Z “h - 

21 315 210 


7512 ® , 1483 ? 9 

155925 1247400 ^ 


23 2 ®° 367 ? _^4 

~102960 ~ 10810800 ^ 

1^82® ? 4 403372 

46305 


2 ®® e 


- / 

403372 ® , 79099 ? g 125621472 ^ 
46305 185220 ^ 91683900 

\ 


4840537 ? 8 
157172400 ^ " 


262 ® In 2 13 ? 0 , 13 2 ®" In 2 13 ? 3 , \ n 

- 2 ® In 2 H-^- 2 ® In 2 

441 441 1323 5670 j 

(-2z‘^ 182981 2 ^ 3753697 ? . 65 2 ® In 2 689 ? 

-- 2 ® 4 -^- 2 ® In 2 

21 92610 5556600 441 13230 

2 2 ® 4 ? 7 13 2 ® 8 ? g 

~ 3465 ~ 10395 ^ 


63 315 ^ 3465 10395 ^ 

712®° ? 44 372®^ ? 43^ 

"^540540 54054 ^ ~ 16216200 ~ 4054050 ^ ) 

/- 22 ® 4 ? _^0 1422 ^ 31 ? 9292 ° 8 ? 

"^ 1^1 135 ^ 51975 ~ 51975 ^ 2702700 96525 ^ 

^5 2 ® 97*^5 9582238912 ® 239560304 ? _^7 

'^ll^'^iiTo^ 6051137400 3781960875 ^ 


6 



Gravitational Waves by a Particle in Cireular Orbits 


11 


1571^6 In z 3142 i 7 


218295 1091475 


z In 2 : e 


-20 2 ^ i 
441 1% 




D7n _ 

-/tc, , - 


2 ^’^ 2 i o 7z^ 
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\^495 1485 19305 15015 
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1621620 737100 ^ 


pin _ 
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2^6 


6 ^, 


.9 2 z 

^ 135135 405405 ^ 


/-32 2 ^ 163 i A 

^I 19305 ~ 405405 ^ 

mn 8 2 ^ 2i 70 

225225 225225 

and the Wronskian for i = 2, and 3 is given by 


oocooc oocooc ’ 


^ 2 " =—1 + 6 


5i TT . , ^ 

y - - - z (7 + ln 2 ) 
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9 
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39 7r (7 + ln 2 ) 5z 2 / 1 r,\ 39 i , , o\2 
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TT (7 + ln2)^ i 3 

+ -^ + -(7 + ln2f 


(A-2) 


where 7 and C(n) represent the Euler constant and the Riemann zeta function, 
respectively. In the expression of we have assumed that u> is positive. The 

expression corresponding to negative lo can be obtained by using the relation = 

For £ > 4, we onR need the forms of valid up to 0(e), which can be found 
in Tagoshi and Sasaki. qI 

Using the above formulas, it is straightforward to obtain the luminosity of the 
gravitational waves. We define 7^3 


as 


dE 1 f dE\ 

^ i,m 


(A-3) 


where {dE/dt)^ = {32/5){fi/M)‘^v^^ is the luminosity given by the well-known 
quadrupole formula. We show r]£^m only for m > 0, since = 'Hi-m ■ 


m,2 = 1 - 


107R 


-h v' 




6 


21 

99210071 
^ 1091475 
19136 vru^ 
1323 


-|- 47ru'^ -|- 


4784 428 vru^ 


1323 


17127 ^ 16 TT^ 
105 3 


21 

3424 In 2 
105 


1712 Inu 
105 


of 27956920577 1831847 1712 366368 In 2 183184 In u' 

+ I_ __ -—— -——- 


81265275 2205 63 2205 

q / 396840284 TT 6848 7 tt 13696 tt In 2 6848 7 r Inu 

-t- (- - - 


2205 


-b V 


10 


-b V 


11 


1091475 
/187037845924 
6257426175 
16380416 In 2 
138915 

-111827682308 71 
81265275 


105 

8190208 7 


-b 


105 

76544712 


105 


138915 3969 

8190208 Inu' 


-b 


138915 

732736771 


1465472 71 In 2 
+-xxxz-+ 


2205 


36 


17 


-b 


71 V 


2205 

732736 71 Inu 
2205 
2215u® 177 iu^ 


(A-4) 


-b V' 


f 15707221 


18 254016 

107 7 Ti^ 


26195400 
2215 71 


945 




252 

107 In 2 
945 


107 Inu 
945 


127008 
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6435768121 , 1819 7 
'57210753600 13230 


17TT^ , 1819 In2 ^ 1819 Ini; 
378 13230 13230 


m,3 = 


715707221TT 21477r 214 tt In2 2147r InuX 

^ V 13097700 945 945 945 ) ’ 

1215 1215 , 3645 TTi;® ^ 243729 v® 3645 

896 112 448 9856 56 


25037019729 473857 , 3645 vr^ 

125565440 1568 224 


47385 In 2 


47385 In 3 47385 lni;\ 

1568 1568 ) 

7311877ri;9 

4928 

10 / 2074855555161 47385 7 3645 vr^ 


- 


47385 In 2 


1381219840 


47385 In 3 47385 \nv\ 

^ 196 196 ) 

/ 75111059187 TT 1421557 TT 


%,2 = 


V 62782720 784 

142155 TT In3 142155 tt In?; \ 

784 784 ) ’ 

5i;^ 193 w® , 207ri;^ ^ 86111?;® 

"63 5^ 63 280665 


142155 TT In 2 
784 


960188809 

178783605 


10407 , 80 vr^ 
1323 189 


772 TT 
' 567 

2080 In 2 


1040 Ini; 


344444 VTi;!! 


280665 


V3,i = 


i;® - 


v‘^ ^ Trv° ^ 437?;® 

~ 1512 4032 266112 

1137077 137 

' 50854003200 ~ 42336 


13 In 2 
42336 


13 In?; 
42336 


437 Tri;^ 


133056 


38943317051 

5034546316800 


-1137077 TT 
25427001600 


137 7r 
21168 


TT^ 13 In2 ^ 13 In?; 
Tim 7938 7938 


137r ln2 13 7 r In?; 


21168 


21168 


^4,4 = 


1280?;^ 151808?;® ^ 10240vr?;'^ ^ 560069632?;® 1214464 tt; 

567 6237 567 6243237 6237 

10 /36825600631808 257392647 8192025739264 In! 

^ I 88497884475 392931 1701 130977 
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25739264 lnv\ 4480557056 -k 
392931 ) ^ 6243237 ’ 

729 28431 v® , 2187 Tri;^ ^ 620077923 85293 

4480 24640 2240 246646400 12320 ’ 

437^6 , 20 7ri;7 ^ 7199152 u® 1748 Tri;^ 

3^ ~ 43659 3969 218513295 43659 

10 / 9729776708 251367 80 vr^ 50272 In 2 25136 In v 

^ I 619485191325 ~ 2750517 11907 ~ 2750517 2750517 


(A-9) 

(A-10) 


+ 

^ 4,1 = 

%,5 = 

%,4 = 
%,3 = 
%,2 = 
%,1 = 

^? 6,6 = 
%,5 = 
%,4 = 
%,3 = 


28796608 
218513295 ’ 

101 ?;® , , 7478267 101 

282240 ~ 4656960 141120 139848508800 ~ 2328480 ’ 

9765625 v® 2568359375w® , 48828125 
2433024 47443968 1216512 

, 706047851562512841796875 tt 
25904406528 23721984 ’ 

4096 w® 18231296 ^ 32768 

13365 6081075 13365 ’ 

2187 u® 1509031;® ^ 6561 ttv® ^ 600654447 452709 

450560 ~ 2928640 225280 2665062400 1464320 

4?;® 15644 u^® ^ 167ri;^^ 

40095 ~ 18243225 40095 ’ 

ti® 179 W® ^ TTV® 

127733760 ~ 2490808320 63866880 

290803 v^® 1797ri;ii 

"^971415244800 ~ 1245404160’ 

26244 u® 2965572 i;i® ^ 314928 
3575 25025 3575 ’ 

244140625?;^® 

435891456 ’ 

131072 t;® 4063232 vi® ^ 1048576 

9555975 22297275 9555975 ’ 

59049 ?;^® 

98658560 ’ 


5733585 


4z;^® ^ 

495495 5733585’ 


%,i = 


m,7 = 


m,5 = 


7192209024’ 
96889010407 ?;i® 
7116595200 
6103515625^1® 
181330845696 ’ 


(A-ll) 

(A-12) 

(A-13) 

(A-14) 

(A-15) 

(A-16) 

(A-17) 
(A-18) 
(A-19) 
(A-20) 
(A-21) 
(A-22) 
(A-23) 
(A-24) 
(A-25) 



Gravitational Waves by a Particle in Cireular Orbits 


15 


1594323 
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